Accepted to ApJ 

Preprint typeset using I^T^jX style cmulatcapj v. 6/22/04 



STELLAR FEEDBACK IN MOLECULAR CLOUDS AND ITS INFLUENCE ON THE MASS FUNCTION OF 

YOUNG STAR CLUSTERS 

S. Michael Fall 1 , Mark R. Krumholz 2 , and Christopher D. Matzner' 5 

Accepted to ApJ 



o 

(N 



0< 

6 

(N 
> 
00 

m 

(N 
(N 

d 
o 



X 



ABSTRACT 

We investigate how the removal of interstellar material by stellar feedback limits the efficiency 
of star formation in molecular clouds and how this determines the shape of the mass function of 
young star clusters. In particular, we derive relations between the power-law exponents of the mass 
functions of the clouds and clusters in the limiting regimes in which the feedback is energy-driven and 
momentum-driven, corresponding to minimum and maximum radiative losses, and likely to bracket all 
realistic cases. We find good agreement between the predicted and observed exponents, especially for 
momentum-driven feedback, provided the protoclusters have roughly constant mean surface density, 
as indicated by observations of the star-forming clumps within molecular clouds. We also consider a 
variety of specific feedback mechanisms, concluding that H n regions inflated by radiation pressure 
predominate in massive protoclusters, a momentum-limited process when photons can escape after 
only a few interactions with dust grains. We show in this case that the star formation efficiency 
depends on the masses and sizes of the protoclusters only through their mean surface density, thus 
ensuring consistency between the observed exponents of the mass functions of the clouds and clusters. 
Our numerical estimate of this efficiency is also consistent with observations. 

Subject headings: galaxies: star clusters — HII regions — ISM: bubbles — radiative transfer — stars: 
formation — stars: winds, outflows 



1. INTRODUCTION 

Most stars form in protoclusters in dense molecular 
clumps (Lada & Lada 2003; McKee & Ostriker 2007). 
The energy and momentum injected by young stars then 
removes the remaining interstellar material (ISM), thus 
ending further star formation and reducing the gravita- 
tional binding energy of the protoclusters. This feedback 
limits the efficiency of star formation — the ratio of final 
stellar mass to initial interstellar mass — to only 20 — 30%, 
and leaves many protoclusters unbound, with their con- 
stituent stars free to disperse. Even those protoclusters 
that survive will lose some stars by ISM removal and 
subsequent processes. 

Two of the best probes of these formation and disrup- 
tion processes are the mass functions of molecular clouds 
and young star clusters, defined as the number of objects 
per unit mass, tjj(M) = dN/dM. For molecular clouds, 
the best-studied galaxies are the Milky Way and the 
Large Magellanic Cloud (LMC), while for star clusters, 
they are the Antennae and the LMC. In these and other 
cases, the observed mass functions can be represented 
by power laws, ip(M) oc M@ , from I0 4 M Q or below to 
10 6 M Q or above. Giant molecular clouds (GMCs) iden- 
tified in CO surveys have [3 w —1.7 (Rosolowsky 2005; 
Blitz et al. 2007; Fukui et al. 2008). This exponent is also 
found for massive self-gravitating clumps within GMCs, 
the formation sites of star clusters, whether they are 
identified by CO emission (Bertoldi & McKee 1992) or 
higher-density tracers such as C 18 0, 13 CO, and thermal 
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dust emission (Reid & Wilson 2006; Munoz et al. 2007; 
Wong et al. 2008). Young star clusters have j3 w —2.0 
(Elmegreen & Efremov 1997; McKee & Williams 1997; 
Zhang & Fall 1999; Dowell et al. 2008; Fall et al. 2009; 
Chandar et al. 2010). The similar exponents for clouds 
and clusters indicate that the efficiency of star formation 
and probability of disruption are at most weak functions 
of mass. This conclusion is reinforced by the fact that 
(3 is the same for 10 7 — 10 8 yr-old clusters as it is for 
10 6 - 10 7 yr-old clusters (Zhang & Fall 1999; Fall et al. 
2009; Chandar et al. 2010). 

These empirical results may at first seem puzzling. 
Low-mass protoclusters have lower binding energy per 
unit mass and should therefore be easier to disrupt than 
high-mass protoclusters. Indeed, several authors have 
proposed that feedback would cause a bend in the mass 
function of young clusters at M ~ 10 5 M Q , motivated in 
part by the well-known turnover in the mass function of 
old globular clusters (Kroupa & Boily 2002; Baumgardt 
et al. 2008; Parmentier et al. 2008). For young clusters, 
such a feature is not observed (as noted above), while 
for globular clusters, it arises from almost any initial 
conditions as a consequence of stellar escape driven by 
two-body relaxation over ~ 10 10 yr (Fall & Zhang 2001; 
McLaughlin & Fall 2008, and references therein). Nev- 
ertheless, we arc left with an important question: What 
are the physical reasons for the observed similarity of 
the mass functions of molecular clouds and young star 
clusters? 

The goal of this Letter is to answer this question. In 
Section 2, we derive some general relations between the 
mass functions of clouds and clusters. In Section 3, we 
review a variety of specific feedback processes and esti- 
mate the star formation efficiency for radiation pressure, 
the dominant process in massive, compact protoclusters. 
We summarize in Section 4. 
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2. MASS FUNCTIONS 

The radiative losses inside protoclusters determine how 
much of the energy input by stellar feedback is avail- 
able for ISM removal. This in turn depends on the cloud 
structure and the specific feedback mechanisms involved, 
but two limiting regimes bracket all realistic situations: 
energy-driven, with no radiative losses, and momentum- 
driven, with maximum radiative losses. We estimate the 
mass of stars M* and the corresponding efficiency of star 
formation, £ — M*/M, needed to remove the ISM from 
protoclusters in these regimes as follows. We character- 
ize a protocluster by its mass M, half- mass radius Rh, 
mean surface density E, velocity dispersion V m (includ- 
ing the orbital motions of the stars and the turbulent 
and thermal motions of the interstellar particles), RMS 
escape velocity V e , and crossing time r c . For simplic- 
ity, we neglect rotation, magnetic support, and external 
pressure (but see Section 3). Then the properties of a 
protocluster are related by V% = 0AGM/R h , V e = 2V m , 
t c = Rh/Vm (Spitzer 1987) and E « (M/2)/{itRl). We 
also assume that the sizes and masses of protoclusters 
are correlated, with a power-law trend, Rh oc M a . 

In Figure 1, we plot E and Rh against M for star- 
forming molecular clumps in the Milky Way, based on 
measurements of CS, C 17 0, and 1.2 mm dust emission in 
three independent surveys (Shirley et al. 2003; Faundez 
et al. 2004; Fontani et al. 2005). These clumps were 
selected for their star-formation activity (water masers, 
IRAS colors), not their surface density. Evidently, there 
is a strong correlation between Rh and M, and almost 
none between E and M, corresponding to a ~ 1/2. The 
typical surface density is close to the value E ~ 1 g cm~ 2 
expected from theory (McKee & Tan 2003; Krumholz 
et al. 2007; Krumholz & McKee 2008). 4 We assume that 
the Milky Way relations also hold in other galaxies and 
extend up to ~ 10 6 M Q , although it is conceivable that 
they break down above ~ 10 5 M Q . Indeed, Baumgardt 
et al. (2008) and Parmentier et al. (2008) assume that 
Rh is not correlated with M (corresponding to a = 0), 
based on observations of gas-free clusters (e.g. Murray 
2009). However, since ISM removal necessarily occurs 
during the earlier, gas-dominated phase, awl/2 seems 
more appropriate in the present context. As we show 
here, a « 1/2 is also needed to reconcile the observed 
mass functions of molecular clouds and star clusters. 

The rates of energy and momentum input are propor- 
tional to the stellar mass 5 : E oc EM and P oc EM. 
We assume that the timescale for ISM removal is a few 
crossing times: At ~ (1 - 10) x t c (Elmegreen 2000, 2007; 
Hartmann et al. 2001; Tan et al. 2006; Krumholz & Tan 
2007). Thus, the total energy and momentum input are 
E EAt oc EMR h /V m and P « PAt oc EMR h /V m . 
These reach the critical values needed to remove the ISM, 
Bait = kMV e 2 and P crit = MV e , for 



£ oc V^/Rh oc m (3 ~ 5q)/2 (energy - driven), (la) 
EctVe /Rh oc M x ~ 2a (momentum - driven), (lb) 

4 For reference, the Larson (1981) relation for CO-selected clouds 
corresponds to a much lower surface density, S ~ 0.02 g cm -2 . 

5 This is a good approximation for all feedback mechanisms ex- 
cept protostellar outflows, which inject energy and momentum in 
proportion to the star formation rate. Outflows, however, are non- 
dominant in massive protoclusters; see Table 1. 
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Fig. 1. — Surface density S and radius R plotted against mass M 
for star-forming molecular clumps from measurements by Shirley 
et al. (2003; circles, CS emission), Faundez et al. (2004; triangles, 
dust emission), and Fontani et al. (2005; squares, C 17 and dust 
emission). We exclude clouds with M < IOOMq, since they cannot 
form clusters. The lines are least-squares regressions (log R against 
log Af) with a = 0.5 fixed (solid) and a = 0.38 ± 0.023 (dashed). 
The true uncertainty on a is undoubtedly larger than the quoted 
one-sigma error. 



For a = 1/2, the efficiency has little or no dependence 
on mass: £ oc M 1 / 4 in the energy-driven regime, £ = 
constant in the momentum-driven regime. For a = 0, 
the variation is much stronger: £ oc M 3 / 2 and £ oc M, 
respectively. These relations are valid for £ < 0.5. 

Any dependence of £ on M will cause the mass func- 
tions of star clusters ip*(M*) and molecular clouds ip(M) 
to have different shapes. For the moment, we confine 
our attention to clusters young enough to be easily rec- 
ognizable even if they are unbound and dispersing. This 
limit is ~ 10 7 yr for extragalactic clusters such as those 
in the Antennae (Fall et al. 2005). In this case, the 
mass functions of the clusters and clouds are related by 
ip*(M*)dM* oc ip(M)dM (with a coefficient greater than 
unity if several clusters form within each cloud). For 
ip(M) oc M and £ oc Af 7 , we have ij)*{M*) oc Aff* with 
/3* = (/3 — 7)/(l + 7). Equations (la) and (lb) then 
imply 



ft. 
ft 



2/3 + 5a - 3 

5(1 - a) 
P + 2a - 1 
2(1 -a) 



(energy — driven), (2a) 
(momentum — driven). (2b) 



These expressions give ft = (3 for a = 3/5 and 1/2, 
respectively. Thus, the similarity of the mass functions 
of clusters and clouds (ft « ft) requires that the latter 
have approximately constant mean surface density (0.5 < 
a < 0.6), no matter what type of feedback is involved. 

Before proceeding, we make a small correction. 
For clouds, the observed mass function ip a (M) repre- 
sents the true mass function at formation ip(M) (i.e., 
the birthrate) weighted by the lifetime: ip a (M) oc 
ip(M)ri(M). We assume, as before, that lifetime is pro- 
portional to crossing time: t; oc r c oc Af( 3Q-1 )/ 2 . Then 
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the exponents of the true and observed mass functions 
are related by /3 — f3 — (3a — l)/2. Inserting this into 
Equations (2a) and (2b), we obtain 



2 (fl, + a- l) 

5(1 - a) 
2/3 Q + a - 1 
4(1 - a) 



(energy — driven), (3a) 
(momentum — driven). (3b) 



We now evaluate Equations (3a) and (3b) with (5 = 
— 1.7, the observed exponent of the mass function of 
molecular clouds (Rosolowsky 2005; Reid & Wilson 2006; 
Munoz et al. 2007; Wong et al. 2008; Fukui et al. 2008). 
For constant mean surface density (a = 1/2), we find 
/?* = —1.8 in the energy-driven regime and ft* = —2.0 in 
the momentum-driven regime. These predictions agree 
nicely with the observed exponents of the mass functions 
of young star clusters, /?* « —2.0 (with typical uncer- 
tainty A/3* w 0.2). Our model is clearly idealized, but 
the scalings, and thus the agreement between the pre- 
dicted and observed /3*, should be robust. For constant 
size (a — 0), however, we find /3* = —1.1 in both the 
energy-driven and momentum-driven regimes, in definite 
conflict with observations. 

The mass function of star clusters older than ~ 10 7 yr 
depends on the proportion that remain gravitationally 
bound. This in turn depends on the efficiency of star 
formation £ and the timescale for ISM removal At rela- 
tive to the crossing time r c . Both analytical arguments 
and TV-body simulations indicate that young clusters lose 
most of their stars for £ < 0.3 and At <C t c but re- 
tain most of them for £ > 0.5 or At > t c (Hills 1980; 
Kroupa et al. 2001; Kroupa & Boily 2002; Baumgardt 
& Kroupa 2007). Thus, as long as £ and At/r c are, on 
average, independent of M, as they are for protoclus- 
ters with constant mean surface density (a = 1/2) and 
momentum-driven feedback, ISM removal will not alter 
the shape of the mass function (although its amplitude 
will decline). This is consistent with the observed expo- 
nents /?* w —2.0 for clusters both younger and older than 
10 7 yr in the Antennae and LMC (Zhang & Fall 1999; 
Fall et al. 2009; Chandar et al. 2010). 

In all other cases, £ increases with M, and a higher 
proportion of low-mass clusters is disrupted, causing a 
flattening or a bend at £ w 0.3 — 0.5 in the mass func- 
tion. The exact shape depends on At/r c , dumpiness 
within protoclusters, and other uncertain factors. If the 
efficiency has a weak dependence on mass, as it does for 
constant mean surface density (a = 1/2) and energy- 
driven feedback (£ oc M 1 / 4 ), the predicted /3* might be 
marginally consistent with observations over a limited 
range of masses (10 4 M Q < M < 1O 6 M ). However, 
for constant size (a = 0), the variations are so strong 
(£ cx M z / 2 and £ oc M) that we expect major differences 
between the mass functions of clusters younger and older 
than 10 7 yr, in clear contradiction with observations. 

Our simple analytical model agrees, at least qualita- 
tively, with the numerical calculations by Baumgardt 
et al. (2008) and Parmentier et al. (2008). They present 
results for energy-driven feedback by supernovae in pro- 
toclusters with uncorrelated sizes and masses. In some 
cases, they find a bend in the mass function of young clus- 
ters at M ~ 1O 5 M , while in others, they find a flattened 
power law with /?* « — 1 (see Figure 4 of Baumgardt 
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Fig. 2. — Feedback in protoclusters of mean surface density 
S and mass M. Radiation pressure is the dominant mechanism 
throughout the shaded region. The lines show where each mecha- 
nism alone achieves E = 0.5. These allow for partial sampling of 
the stellar IMF and hence differ slightly from the power laws in 
Table 1 (noticeable only for M < KTMq). 



et al. 2008). As we have already noted, these results are 
expected for a = 0, and they are inconsistent with the 
observed mass functions of young clusters. 

3. STAR FORMATION EFFICIENCY 

We now consider five specific feedback mechanisms: 
supernovae, main-sequence winds, protostellar outflows, 
photoionized gas, and radiation pressure. For the first 
four, we review results from the literature. Supernova 
feedback begins only after the > 3.6 Myr lifetimes of 
massive stars. Unless turbulence within a protoclus- 
ter is maintained by feedback or external forcing, stars 
would form rapidly and consume its ISM, with £ — > 1 in 
1 — 2 crossing times. This implies that supernovae can 
dominate only for 2r c > 3.6 Myr unless another mecha- 
nism somehow keeps £ small without expelling much ISM 
(Krumholz & Matzner 2009). However, even in this con- 
trived situation, supernovae would play only a secondary 
role. Main-sequence winds are not effective if their en- 
ergy is able to leak out of the bubbles they blow (Harper- 
Clark & Murray 2009). As a result of this leakage, winds 
simply provide an order-unity enhancement to radiation 
pressure (Krumholz & Matzner 2009). Protostellar out- 
flows can only remove the ISM from protoclusters with 
escape velocities below about 7 kms -1 (Matzner & Mc- 
Kee 2000). Photoionized gas is important as a feedback 
mechanism only when its pressure exceeds that of radi- 
ation throughout most of an H n region. This in turn 
requires that the H II region be larger than the radius 
r ch at which — -Pgasj a condition harder to satisfy 
in massive, compact protoclusters (Krumholz & Matzner 
2009). 

We summarize these results in Table 1 and Figure 2. 
As the plot shows, the mechanisms discussed thus far 
are relatively ineffective in protoclusters with M > 10 4 
M and E > 0.1 g cm~ 2 . We therefore turn to radi- 
ation pressure. This would be an energy-driven feed- 
back mechanism if all photons, even those re-radiated 
by dust grains, remained trapped within a protoclustcr. 
However, this is possible only if the protoclustcr is so 
dense and smooth that the covering fraction seen from 
its center exceeds ~ 90 % in the infrared (Krumholz & 
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TABLE 1 

Feedback Mechanisms 



Mechanism 



Type 



Limitation 



Threshold* 



Evaluated* 



Supernovae 
Main-sequence winds 
Protostellar outflows 
Photoionizcd gas 
Radiation pressure 



Energy 
Either 3 - 
Momentum 
Momentum 
Momentum 



Too late 
Relatively weak a 
Confined in massive clusters' 3 
Crushed by P rad c 



t c Ri 1.8 Myr 

Never 
V e « 7 km s- 1 
5 49 ~ 21R h /pc 
Equations (6) and (7) 



ll' 3 

> 0.17M" 1 
S « O.ISM^ 1 
So « 1.2 



So 



^Parameters required for £ — 0.5. Evaluations assume a fully-sampled stellar IMF. Notation: S49 = 5/10 49 s _1 
(ionization rate), M 4 = M/10 4 M©, S = S/g cm -2 . 

a Stcllar winds are energy-driven and dominant if trapped, but are expected to leak, making them momentum-driven 
and weak. 

b Based on Equation (55) of Matzncr & McKec (2000), updated with f w v w = 80 km s" 1 (Matzncr 2007). 

c Based on Equation (4) of Krumholz & Matzncr (2009) for the blister case, with the coefficient reduced by a factor of 
2.2 to correct an error in the published paper and updated with (L/M.) = IUOLqM' 1 and (S/M,) = 6.3 x 10 46 s" 1 
M" 1 (Murray & Rahman 2010). 



Matzner 2009). More realistically, the protocluster would 
be porous enough that photons could escape after only a 
few interactions with dust grains, and radiation pressure 
would then be a momentum-driven feedback mechanism. 
The following analysis extends that of Elmegreen (1983), 
Scoville et al. (2001), Thompson et al. (2005), Krumholz 
& Matzncr (2009), and Murray et al. (2010). 

We consider an idealized, spherical cloud of mass M 
and outer radius R, with an internal density profile 
p oc r~ k (hence R h = 2~ 1 /^~ k ^ R). Radiation from 
young stars near the center of the cloud ionizes the gas 
and drives the expanding outer shell of this H n region. 
After a time t, the momentum imparted to the shell is 
Ps = ftr&pLt/c, where L is the stellar luminosity (as- 
sumed constant for simplicity), and / trap ~ 2 — 5 ac- 
counts for assistance from main-sequence winds and in- 
complete leakage of starlight and wind energy (Krumholz 
& Matzner 2009). Neglecting gravity for the moment, the 
velocity and radius of the shell are related by v s = r\r s jt 
with r] = 2/(4 — k). Thus, when the shell reaches the 
cloud surface (r s = R), it has swept up all the remaining 
ISM, with mass M, 

by 



g — (1—£)M, and has a velocity given 

VftrapLR 



(4) 



c(l-£)M' 

We specify the condition for ISM removal against grav- 
ity of the protocluster by v 2 (R) = a cr i t GM/(5R), where 
a cr it is a parameter of order unity that accounts for mag- 
netic support and other uncertain factors (discussed be- 
low). The required luminosity, from equation (4), is 

a clit Gc(l - £)M 2 



L = 



(5) 



5?7/tra pj R 2 

The fundamental scaling L oc (M/R) 2 oc arises here 
in the same way it does for the growth of supermassive 
black holes and galactic spheroids (Fabian 1999; King 
2003; Murray et al. 2005). Rewriting Equation (5) in 
terms of X = M/(irR 2 ) and M* = £M and solving for 
£, we obtain our basic result 

£ = ^^, (6) 



with 



577/trap(£/M.) 

na crit Gc 



1.2 



/trap \ 
«crit / 




gem 



(7) 



Fig. 3. — Star formation efficiency £ as a function of mean surface 
density E, computed from Equations (6) and (7) with the indicated 
values of /t r a P /a C rit- 



The coefficient in the last equation is based on 77 = 2/3 
and L/M* = 114OL /M (see notes to Table 1). Re- 
gardless of the exact value of /trap/«crit, we note that 
£ depends on M and R only through S. Thus, when 

5 is constant, £ is independent of M, and the mass 
functions of clusters and clouds have the same exponent 

(/?* = P « Pa). 

Figure 3 shows £(E) computed from Equations (6) and 
(7). Clearly, £ increases monotonically with E from to 
1, reaching £ = 0.3 for E - 0.5(/ tra p/acrit) g cm" 2 . 
We expect / trap <~ a cr i t ~ 2-5. The escape velocity 
from the surface of an unmagnetized cloud corresponds 
to a cr ;t = 10, while the internal velocity dispersion, pos- 
sibly sufficient for some ISM removal, corresponds to 
a cr ; t » 1,3. A shell driven by a constant force requires 
a cr it = 2.3 (for k = 1; see Equation (A17) of Matzner 

6 McKee (2000)). We consider a crit « 2 to be plau- 
sible; certainly a protocluster boils violently and loses 
mass rapidly using this condition. Our intent here is not 
to make a detailed comparison between the model and 
observations. Given the simplicity of the former and the 
uncertainties in the latter, it is gratifying that they agree 
even roughly with each other. 

4. CONCLUSIONS 
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This Letter contains two main results. The first is the 
relation between the power-law exponents of the mass 
functions of molecular clouds and young star clusters, 
fi a and /?*, in the limiting regimes in which stellar feed- 
back is energy-driven and momentum-driven, Equations 
(3a) and (3b), which bracket all realistic cases. The pre- 
dicted /?* depends significantly on the initial size-mass 
relation of the protoclusters. We find good agreement 
between the predicted and observed /?*, especially for 
momentum-driven feedback, for S oc M/Rf « constant, 
the relation indicated by observations of gas-dominated 
protoclusters. In this case, the star formation efficiency 
is independent of protocluster mass, ensuring that the 
fraction of clusters that remain gravitationally bound fol- 
lowing ISM removal is also independent of mass. 

The second main result is an estimate of the star for- 
mation efficiency in protoclusters regulated by radiation 



pressure, Equations (6) and (7). This is likely to be 
the dominant feedback process in massive protoclusters. 
We show that £ depends on M and Rh only through 
the mean surface density S, which in turn guarantees 
consistency between the observed power-law exponents 
of the mass functions of molecular clouds and young 
star clusters according to our general relations. For 
S ~ 1 g cm~ 2 , we estimate £ ~ 0.3, in satisfactory 
agreement with observations. 
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